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Abstract 

We model the dynamics of Magellanic Stream with the ram-pressure scenario in 
the logarithmic and power-law galactic halo models and construct numerically the 
past orbital history of Magellanic clouds and Magellanic Stream. The parameters 
of models include the asymptotic rotation velocity of spiral arms, halo flattening, 
core radius and rising or falling parameter of rotation curve. We obtain the best fit 
parameters of galactic models through the maximum likelihood analysis, comparing 
the high resolution radial velocity data of HI in Magellanic Stream with that of 
theoretical models. The initial condition of the Magellanic Clouds is taken from 
the six different values reported in the literature. We find that oblate and nearly 
spherical shape halos provide a better fit to the observation than the prolate halos. 
This conclusion is almost independent of choosing the initial conditions and is valid 
for both logarithmic and power-law models. 

Key words: galaxies: Magellanic Clouds, Galaxy: structure 
PACS: 98.56.Tj 



1 Introduction 



Analysis of rotation curve of spiral galaxies reveals a faster rotation curve 
at the outer parts of spiral arms compare to what is expected based on the 
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distribution of visible matter. Hence invoking a halo shape dark component for 
galaxies is essential to explain the dynamics of galaxies (Bosma 1981; Rubin 
and Burstein 1985). Despite using a large number of techniques, the shape of 
the galactic halo is not well known and is still in debate. It is even not clear 
if the halo is prolate or oblate. Most halo models expect significant deviations 
from the spherical symmetry (Frenk et al. 1988; Dubinski and Carlberg 1991; 
Oiling and Merrifield 2000 and Ruzicka et al. 2007). An accurate and reliable 
determination of the shape and mass distribution of the galactic halo requires 
the analysis of the kinematic properties of several halo objects (i.e. satellite 
galaxies and globular clusters). Another promising method is studying the 
stellar and gaseous streams which produced by disruption and accretion of 
small galaxies orbiting around the Milky Way (MW). The shape and the 
kinematics of such streams are strongly influenced by the overall properties 
of the underlying potential. There are several examples of such streams which 
have been discovered in recent years (Odenkirchen et al. 2001; Grillmair & 
Dionatos 2006; Ibata et al. 2001; Newberg et al. 2002; Majewski et al. 2003). 
The most earliest discovered stream is the Magellanic Stream (MS) which is 
a gaseous streams associated to the Magellanic Clouds (MCs). 



The MS is a narrow band of neutral hydrogen clouds lies along a great circle 
from (I = 91°, b = -40°) to (I = 299°, b = -70°), started from the MCs and 
oriented towards the south galactic pole. The radial yelocity and the column 
density of MS ha s been measured by several groups (IMathewson et al. 19871 ; 



Brians et al. 20051 ). Observations show that the radial velocity of MS with 
respect to the Galactic center changes from kms~ l at the nearest distance 
from the center of MCs to —200 kms -1 at the tail of stream. Another feature 
of this structure is that the column density falls about two orders of magnitude 
along the stream. The two main competing scenario for explaining the origin 
of MS are based on the tidal interaction and the ram pressure stripping which 
predict different ages and spatial distributions for MS. In the tidal model the 
idea is the gravitational stripping of gas from the MCs due to the strong 
tidal interaction of LMC and SMC with Galaxy at their perigalactic passage 
(Murai and Fujimoto 1980; Lin & Lynden-Bell 1982; Gardiner & Noguchi 
1996; Connors et al. 2006). This leads to the formation of tidal tails emanating 
from the opposite side of MCs. In the tidal scenario, although the observed 
radial velocity profile of the stream has been modeled remarkably well, the 
smooth HI column density distribution does not match with the observations 
and the expected stars in the stream have not been observed yet. The other 
hypothesis is the ram pressure stripping of gas from the MCs by an extended 
halo (Moore & Davis 1994; Heller & Rohlfs 1994; Sofue 1994). In this model, 
there is a diffused halo around the Galaxy which produces a drag on the gas 
within the MCs and causes material to escape and form a trailing gaseous 
stream. Compared to the tidal scenario, the ram pressure model allows a 
better production of HI, compatible with the column density of MS. 
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In this work, we study the dynamics of MS for two generic galactic mod- 
els to constraint t heir parameters based on detailed HI observations of MS 



( iBriins et al. 20051 ). The main problem with modeling the interaction of the 
MCs with the MW is the uncertainty in the initial condition of MCs. Here 
we do calculations for six different initial conditions of MCs which have been 
reported in the literatures (Tabled]). 

The paper is organized as follows: In section [2] we introduce the potential 
models of MW and MCs. In section [3] we study the dynamics of MS in various 
Galactic models and compare them with the observations. The best parame- 
ters of model is obtained through the maximum likelihood analysis in section 
HJ The conclusion is given in section [51 



2 Models for Milky Way and Magellanic Clouds 



2.1 Milky Way 



The galactic model for the gravitational potential of an spiral galaxy consists 
of a disk, bulge and halo components. For the halo component of MW we 
take two sets of models so-called "power-law" and "logarithmic" potentials. 
The axisymmetric power-law and logarithmic models for the galactic halo are 
given by the following potentials (jEvans 1993 fc 19941 ): 



- VaHogiR^ + tf + ziq- 



'V 



and 



(R c 2 + R 2 + z 2 q- 2 )P /r 



(2) 



where ^f a — is the central potential, R c is the core radius and q is the flat- 
tening parameter, q — 1 represents a spherical halo and q ^ 1 gives an elliptical 
shape. The parameter f3 determines whether the rotation curve asymptotically 
rises (/3 < 0), falls (/3 > 0) or is flat (/3 = 0). These models can reproduce the 
observed rotation curve of MW. 



The galactic disc and bulge are the luminous components of the MW. For 
simplicity in calculation at the distances r > 50kpc, we take the gravitational 
contribution of bulge and disk are as poin t like objects with the masses of 3.5 x 
1O 1O M and 1.2 x lO n M , respectively flHelmi 20041 : iRead and Moore 20051 ). 
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2.2 Magellanic Clouds 



For LMC and SMC we adopt the Plummer potential as: 



Gm LtS 



= T7 ^T¥2 V (3) 

[(r - r Li5 ) 2 + KiJ 



where tl,s are the position of the Clouds relative to the MW center, and Kl = 
3kpc and Ks = 2kpc are the core radii. We assume a mass of tul = 2 x 1O 1O M 
for the LMC, and m s = 3 x 1O 9 M for the SMC. 



3 Dynamics of the Magellanic system 



The simulation of MCs requirers the knowledge of 3D orbits of LMC and 
SMC. We observe only a 2D projection on the sky at the present time so 
3D distribution of positions and velocities are partly unknown. While the 
spatial location of MCs is well known, the the velocity field corresponds to this 
structure due to the uncertainty in the transverse velocity is less constrained. 
It is therefore very difficult to construct a three-dimensional orbit that matches 
all observations. 

In this section we use the RK4 method"*! to calculate the equations of motion 
of MCs. We integrate the equation of motion numerically and extract the 
dynamics of MCs by backtracking orbits from their current positions and 
velocities to t = —5Gyr. The equations of motion of the Clouds is given by: 

= g^- s i<PsA\rL,s ~ r^l) + MM)] + *l,s + F drag , (4) 



where 4>s,l represents mutual gravitational interaction of LMC and SMC, (f>c 
is the gravitational potential of Galaxy, f is dynamical friction force, r is the 
distance of each clouds from the Galactic center and F^rag is the hydrody- 
namical drag force. At the first step, before studying the dynamics of MS we 
examine the effect of various friction forces on the dynamics of MCs. 



Runge-Kutta 
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3.1 Evaluation of friction forces 



In this part, we evaluate the contribution of friction forces in the dynam- 
ics of MCs. We consider a simple isothermal hal o where the potenti al is 



given by <pG{ r ) — —v c 2 lnr, and v c = 220kms 1 (IConnors et al. 20061 ). In 



this calculation we take the initial conditions of MCs as reported in HR94 



(jHeller and Rohlfs 19941 ). see Table [TJ In this calculation we use the Cartesian 
coordinate system where the x-axis is defined in the direction from the sun 
to Galactic center, the y-axis is points towards the direction of the circular 
motion of Sun around the Galaxy, and the z-axis is directed to the northern 
Galactic pole. 

(a) Dynamical friction: Due to the Galactic extended massive halo, the satellite 
galaxies passing through halo will be slowed down by the gravitational force 
of the halo so-called dynamical friction. This force makes MCs to spiral slowly 
toward the Galactic c enter. We adopt the followin g standard form of dynamical 



friction for the halo (IBinney and Tremaine 19871 ): 



Gm 2 T o Vr ej 

f L , s = 0.4281nA 5 

r 2 v LyS 



where In A ~ 3 is the Coulomb logarithm. To evaluate the strength of the 
dynamical friction force versus the gravitational force of the halo, we compare 
their corresponding time scales. The dynamical time scale can be given by 
the orbital period as T dyn ~ r/v c , where v c is the circular velocity and r 
is the distance of MCs from the MW. On the other hand the time scale of 
dynamical friction is in order of Tf ric ~ r 2 v c /GM MCs , where M M q s is the mass 
of the Magellanic Clouds. The ratio of Td yn /Tf r i c ~ M.MCsjM-MW shows that 
Tfric 3> Tdyn- Figure ([1^) shows the effect of dynamical friction on the orbit of 
MCs where the spiraling of clouds get faster by increasing the mass of clouds. 

(b) Hydrodynamical drag force: The hydrodynamical friction force results from 
the direct collision of the halo gas with the MCs while passing through the 
halo. The presence of hot diffused gas in the Galactic halo in hydrodynamical 
equilibrium with the dark halo has been proposed by White and Frenk (1991). 
To explain some ionization features discovered in MS, the baryonic halo gas 
should have a temperature in or d er of ~ 10 6 K distrib uted at a radius larger 



than 70kpc fjSembach et al 20031 ; iPutman et al. 20031 ). Constraints from the 



dynamics and the thermal observations propose the density of gas to be in the 
range of 10~ 5 to 10~ 4 cm~ 3 . Assuming MCs as a dense sphere moving through 
the halo, the drag force on this sphere is given by: 

Fdrag = \c ' d (R e ) P gVfr D 2 , (6) 
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Fig. 1. The effect of friction forces on the orbit of MCs. (a) The effect of dynamical 
friction for different value of MCs's mass, (b) The effect of hydro dynamical drag 
force on the orbit of MCs for different value of MCs's size, (c) The effect of hydro- 
dynamical drag force on the orbit of MCs for different value of halo density, (d) 
Comparison of orbits for different types of friction forces with that of frictionless 
orbit. Here is shown that the effect of dynamical friction is more significant than the 
hydrodynamical drag force. The duration of evolution is 5Gyr. The mass of MCs is 
in terms of 1O 1O M0. The initial condition selected as HR94. The plus sign denotes 
the Galactic position at (0,0,0). The black spot denotes the present position of MCs. 

where Cd{R e ) is the coefficient of drag force which is a function of Reynolds 
number R e , p g , the density of halo, v c , the relative speed with respect to 
the gas, and D, the size of the Clouds. The time scale of the hydrody- 
namical friction force can be given by Td rag ~ MmcsV c / Fdrag- which im- 
plies T drag /T dyn ~ M MCs / p g rD 2 . Substituting the parameters as r ~ 50kpc, 
D ~ 5kpc, number density of halo ~ 10~ 4 5 atom /cm 3 and Mmcs ~ 1O 1O M , 
we obtain T drag 3> T^ yn . Figures (EDd) and (Jlfc) show the effect of friction forces 
appeared in equation (jlj) on the orbital motion of MCs. Due to the dynamical 
friction, MCs lose a few percent of total energy in each passage. The effect of 
the gaseous halo density and size of clouds on the magnitude of drag force is 
indicated in Figure ([!]). The apogalactic distance of clouds decreases gradually 
during the evolutionary period. Finally, we calculate the dynamics of MCs in- 
cluding all forces, represented in Figure (EDi). Amongst the dissipative forces, 
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Table 1 

Orbital parameters of LMC used in this work. The position is in kpc and velocities 
are in kms^ 1 and is measured in Galactocentric coordinate. References: Gardiner & 



Noguchi (1996, GN96), Heller & Rohlfs (1994, HR94), Kallivayalil et al. (2006a,b, Kl 
and K2), Mastropietro et al. (2005, M05), and van der Marel et al. (2002, vdM02). 

the dynamical friction force is the dominant term and can deviate the orbit 
of MCs from that of frictionless one. 



3.2 Modeling of the Magellanic Stream 



In order to determine the evolution and the formation history of MS, it is 
essential to know the spatial location and velocities of MCs at the present 
time as the initial condition. As we mentioned before this initial condition 
suffers from the uncertainties mainly from the tangent velocity field of this 
structure. Different initial conditions show significantly different orbits and 
location of the interaction between LMC and SMC. In this regards, recently an 
extended analysis of the parameter space for the interaction of the Magellanic 
system with the MW has been done by Ruzicka et al (2007), using the genetic 
search algorithm combined with an approximate restricted N-body simulation 
(instead of fully self-consistent simulation). Also in another orbital analyzing, 
different sets of initial conditions have been applie d to determine the first 
passage of LMC around the MW (IBesla et al. 20071 ). One of the main aims 
of this work is seeking the results that are common for different selection of 
initial conditions. 



We use the six different sets of initial conditions listed in Table ([T]) and for each 
set, we find out the best space parameters of MW halo. Since the obtained 
parameters depends on the formation mechanism of the MS, we compare our 
results with the other studies which used other formation mechanism. For 
modeling the MS, we follow the continues ram pressure stripping mechanism 
introduced by Sofue (1994) including the dynamical and hydrodynamical fric- 
tion forces. In that work it is assumed that if LMC and SMC approach to 
each other, this sever encounter will most likely disrupt the two galaxies to 
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form the Magellanic bridge (a stream of HI gas between the LMC and SMC). 
The Magellanic bridge within SMC and LMC then has been stripped off by 
the ram-pressure when the MCs was moving through a hot MW halo. Ac- 
cording to the observa tions, at the prese nt time, the gaseous stuff is leaving 



the Magellanic bridge (jBriins et al. 20051 ). We assume that the gas is initially 



distributed inside the MCs. The relative motion of MCs inside the Galactic 
halo makes a drag force on the gas of Magellanic bridge and due to the larger 
cross section of gas compare to the stars, the gas exit from the MCs. The 
released gas slows down by the friction force and accrete into the galactic disk 
and makes the MS. The equation of motion of MS particles can be written as 

- = -^[^MCsi^MS - *MCs\) + <I>G(\TMS\)] + f + Fdrag, (7) 



where tms and rues are the distance of MS and MCs from the center of 
Galaxy, f is the dynamical friction force and Fd ra g is the hydrodynamical drag 
force. As a simple model, MS is considered as a series of sph erical clump s with 
the size of about 0.5kpc and the mass of m = 3 x 1O 6 M ( ISofue 19941 ). The 
first part of the MS (i.e. end of tail) is generated at the time of LMC-SMC 
close approach about 0.5Gyr ago and subsequent clumps is released from the 
MCs after this time. Since the dynamical friction force is proportional to the 
mass of the object, for these clumps, we can ignore this term in comparison 
with the drag force. Base on such simple model, we were able to produce a 
large number of trailing tail across the entire range of halo parameters values. 
By comparing the created MS profile with the observed data we find out the 
best set of halo parameters. In the next section we calculate the orbit of MS 
for two main category of the power-law and logarithmic models. 



4 Results and discussion 



In this section we calculate the equations of motion of MCs using the gravita- 
tional potential of disk, bulge and dark halo. We extract the dynamics of MS 
and obtain the trajectory of MS in various Galactic models. Then we repro- 
duce the radial velocity distribution of HI in MS for different initial conditions 
of MCs (reported in Table [TJ to compare with the observed data. For each 
initial condition, the best fit orbital parameters have been reported in Table 
(J2J). We assumed that the LMC and SMC form a binary system that has been 
in a slowly decaying orbit about the MW for roughly a Hubble time. One 
of the crucial points in studying the dynamics of MS is the location of the 
close approach and the corresponding time which depends on the velocities 
and positions of MCs at the present time. For the heavy halos, the collision 
happens at the closer distance to the present position of MCs than the light 
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Table 2 

The best Parameters of logarithmic model for various initial conditions. The first 
column represents the corresponding reference for the initial condition of MCs. The 
parameters in the second and third columns (V c and q ) have been introduced in the 
text. The fourth column, T co \ is the last close approach of LMC and SCM, the fifth 
column d m i n indicates the distance between LMC and SMC at last close approach. 
The last column shows the minimum x 2 corresponds to each initial condition. The 
mean value of V c and q from the various initial conditions for MCs is q ~ 0.9 and 
V c ~ 187. 

mass halos. 

In what follows we use generic power-law and logarithmic halo models and 
extract the best values for the parameters of the model through the maximum 
likelihood analysis. We use x 2 fit to compare the deviation of the observed 
radial velocity from the theoretical predictions. 

4-1 Logarithmic model 

We use radial velocity of MS to compare the theoretical model with the ob- 
served data. The radial velocity of MS is calculated by averaging over the data 
at each point of MS along the line of sight. The corresponding error bar results 
from the velocity dispersion of the structure which results from the stochastic 
motions of gas. Figures ([2]) represents the relative distance of LMC and SMC 
for six different initial conditions and the radial velocity of MS as a function 
of angular location for the best free parameters of the model. It should be 
noted that radial velocity is measured with respect to an observer located at 
the center of Galaxy. The zero angular distance corresponds to the position 
of MSI (the nearest part of MS to MCs). 

For each set of halo parameters, the radial velocity profile of HI clumps in MS is 
numerically calculated and compared with the observed radial velocity profile. 
The analysis repeated for various initial conditions of the clouds. The best 
value of model parameters (q and V c ) for various initial conditions are given 
in Table El Except HR94 which prefers more oblate halo, for the other sets of 
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Fig. 2. Relative distance of LMS and SMC as a function of time (top) and radial 
velocity versus angular distance along the MS (bottom) for the best-fit logarith- 
mic halo potential. Open squares show the observed radial velocity and curves 
correspond to that of best fit with various initial conditions. 

initial conditions, better agreement between the model and observational data 
is achieved for nearly oblate or spherical halo. The mean value of the best fit 
parameters, averaging over the different sets of initial conditions are q ~ 0.9 
with V c ~ 187A:ms _1 . These values are in agreement with the flattening of the 
MW halo potential obtained by Helmi (20 04) and are also compatible with the 



implied value of circular velocity of MW (IMcgaugh 20081 ; IXue et al. 2008 I ). 



Finally we present the projected trajectory of LMC, SMC and MS in Figure (j3J) 
for various initial conditions. The MCs move on rosset-shape orbits except the 
case of HR94 and K2 which have the banana-shape trajectory. The apogalactic 
distance of LMC decreases gradually during the evolutionary period, which 
reflects the effect of dynamical friction. According to Figure ([3]), for two initial 
conditions represented by K2 and Kl-mean, the spatial extension of MS is 
shorter than the observed size, while for the case of HR94 it is longer than 
the observed size. This result arises from the longer time for the last close 
approach (T co i) in HR94. The minimum relative distance of the LMC and SMC 
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Fig. 3. Y — Z plane orbital path of LMC (red solid line), SMC (dashed line) and 
MS (blue solid line) for the best-fit logarithmic halo potential with various initial 
conditions. The orbital path of the LMC is not significantly affected by the SMC, 
however SMC's path is deviated by LMC. MW is located at (0, 0). The black circle 
denotes the position of last close approach of LMC and SMC. The duration of 
orbital evolution is bGyr. 

at the last close approach (d min ) and the corresponding time is calculated and 
presented in Table [2j 



4-2 Power-law model 



Following the logarithmic model we calculate the radial velocity distribution 
of MS in the power-law galactic model and compare it with the observed data. 
In the power-law model we have 4 free parameters: the asymptotic velocity, V a , 
the halo flattening, q, the core radius, R c and the rising or falling parameter 
of rotation curve (3. The best value of parameters for each initial conditions 
are given in Table El The mean value of the best fit parameters averaging over 
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Table 3 



The best parameters of power-law model for each sets of initial condition. The 
parameters of this model is introduced in equation ([_]). The mean value of the best 
fit parameters averaging over the different set of initial conditions are q ~ 0.85 and 
V a ~ 168, R c ~ 17, j3 1.7. 



the different set of initial conditions are q ~ 0.85, V a ~ 168, R c ~ 17, and 
/3 ~ —1.7. Again similar to the logarithmic halo model, a better agreement 
with observational data is generally achieved for nearly oblate or spherical 
(g < 1) halo model. 

Figure (j_]) shows the time variation of the relative distance of LMC and SMC 
for the last 5 Gyr and radial velocity in Galactic frame as a function of angular 
separation of MS with respect to MSI. Different curves corresponds to various 
initial conditions of the MCs. The projected trajectory of LMC, SMC and MS 
for the best values of the parameters of model is shown in Figure (jSJ). Similar 
to the logarithmic model, in the case of HR94, the generated MS extends 
longer than the observation. From the equations of motion, the minimum 
relative distance of the LMC and SMC at the last close approach and the 
corresponding time is also calculated (see Table _0). 

In order to visualize the compatibility of simulated MS with the observational 
distribution of the stream, we plot the projected spatial distribution of MS 
corresponding to the different initial conditions in Figure ((_]). We find little 
deviation in the projected distribution of MS with the theoretical predictions. 
However, there is a significant difference between the observation and the 
theoretical prediction in the case of HR94 and M05 for both halo models. 
Although the GN96 initial condition could fit with the MS in logarithmic 
model, it no longer traces the MS in power-law model. Figure (EI) shows that 
in the logarithmic model the projected MS using the GN96, vdm02, Klmean 
and K2 initial conditions are nearly indistinguishable and almost overlaying 
HI data of MS. The similar results valid for power-law model except for the 
case of GN96. 
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Fig. 4. Relative distance of LMS and SMC as a function of time (top) and radial 
velocity versus angular distance along the MS (bottom) for the best-fit Power-law 
halo potential. Open squares show the observed radial velocity and curves corre- 
spond to that from the best fit with various initial conditions. 

5 Conclusion 



Summarizing this work, we studied the formation of Magellanic Stream as 
a consequence of interaction of the Magellanic Clouds and MW. We used a 
continuous ram-pressure model for simulating MS. The dynamics of MS is used 
to put constrain on the shape of Galactic halo. In this model the particles of 
halo sweep the hydrogen gas of the Magellanic Clouds while moving through 
the halo. At the first step we reconstruct the orbital motion of the Magellanic 
Clouds as a three body system of Galaxy, LMC and SMC taking into account 
the initial condition of last two objects at the present time. In this scenario, 
at the last close approach of LMC and SMC a part of their gas has been 
ejected from the system to form Magellanic stream. While the Magellanic 
Clouds gradually spiral towards the Galaxy, the stream of gas releases from 
this system and makes Magellanic Stream. We showed that the drag force of 
halo on the stream is larger compare to that on the Magellanic Clouds while 
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Fig. 5. The Y — Z plane orbital path of LMC (red solid line), SMC (dashed line) 
and MS (blue solid line) for the best-fit power-law halo potential with various initial 
conditions. The orbital path of the LMC is not significantly affected by the SMC, 
however SMC's path is deviated by LMC. MW is located at (0, 0). The black circle 
denotes the position of last close approach of LMC and SMC. The duration of 
orbital evolution is bGyr. 



the dynamical friction force on Magellanic Clouds is larger than on Magellanic 
Stream. These friction forces causes the orbit of these structures to spiral faster 
toward the Galaxy. Using this simple model for generating MS, we were able 
to find the best parameters of halo model, comparing the radial velocity of 
Stream with the theoretical model. 



In order to see the effect of different shape of Galactic potentials, we took 
two generic logarithmic and the power-law potentials for the Galactic halo 
and calculate the dynamics of Magellanic Stream for these potentials. Since 
there is a large uncertainty in the initial conditions of MCs at the present 
time, we applied six different initial conditions reported in the literature 
(jBesla et al. 20071 ). Finally we compared the numerical results of radial ve- 
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Fig. 6. HI distribution of MS from the data ( Briins et al. 20051 ) is plotted as a 
polar projection for both logarithmic and power law model. The axis are the galactic 
longitude and latitude as defined in the axis label. The blue dot indicates the location 
of the LMC. The colored lines are the projection of spatial distribution of the MS 
which theoretically determined by several initial conditions under logarithmic and 
power-law MW model. Note that the different initial conditions almost does match 
the data. The HR94 result is indistinguishable from that of the M05 result and 
both deviate from the location of the tip of MS. The similar situation occur for the 
vdm02, Klmean and K2 velocities. 



locity profile of MS with the observed data. Using the likelihood analysis we 
found that Galactic halo is nearly oblate almost in all the initial conditions 
for the two galactic halo models. Furthermore, the preferred value for the 
flattening parameter in our analysis is in agreement with the recent analysis 
which has emp loyed a tidal model fo r the MS and logarithmic potential for the 
Galactic halo (jRuzicka et al. 20071 ). In addition, the corresponding value for 
circul ar velocity is com patible with the recent measurement based on SDSS 
data (IXue et al. 2008 I) and is in agreement with the rotation curve analysis 
flMcgaugh 20081 ). 



We would like to thank C. Briins provided us recent data of Parkes HI survey 
of Magellanic System and his useful comments. 
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